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Abstract. Employing an effective field theory approach to inflationary perturba- 
tions, we analyze in detail the effect of curvature-generated Lagrangian operators on 
various observables, focusing on their running with scales. At quadratic order, we 
solve the equation of motion at next-to-leading leading order in a generalized slow-roll 
approximation for a very general theory of single-field inflation. We derive the result- 
ing power spectrum, its tilt and running. We then focus on the contribution to the 
primordial non-Gaussianity amplitude /nl sourced by a specific interaction term. We 
show that the running of Jnl can be substantially larger than what dictated by the 
slow-roll parameters. 
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1. Introduction 



Inflation [H [2] stands as one of the main pillars of modern cosmology as it naturally 
solves the so-called flatness, horizon and monopoles problems. It also explain the pro- 
duction of density perturbations in the early Universe which then lead to Large Scale 
Structures [3]- [7] in the distribution of galaxies and temperature anisotropies in the Cos- 
mic Microwave Background (CMB) [8]-[l2]. In recent years, many different inflationary 
models have been studied and the advent of new generation experiments (the launch 
of the Planck satellite [131 E] an d the continued analysis of WMAP data [15]) pushes 
further the research in this direction. Indeed, one can now hope that with improved 
sensitivity to deviations from Gaussian statistics, we might be able to probe deeper into 
the dynamics of the inflaton [16]. This translates into including and studying 3-rd and 
higher order (self) interactions in the inflationary Lagrangian [TTJ 13 [191 120] . A detailed 
analysis has been done for a huge variety of models (see [21] for a comprehensive and 
updated review) with precise predictions for the amplitude and the shapes of correlation 
functions. In fact, in performing these investigations, one is aiming at predictions for 
observables such as the scalar spectral index, its running, the primordial bispectrum 
amplitude, Jnl, and its running. The work we present here fits in this context and its 
goal is to highlight the contribution to these observables coming from specific opera- 
tors of a very general effective field theory describing inflation driven by a single scalar 
degree of freedom. A formalism that allows a unified approach to inflation, by means 
of effective field theory techniques, has been presented in [221 E3]. The prescription 
of [23] is appealing for several reasons: one can see there is a clear cut dictionary be- 
tween a given inflationary model and a specific linear combination of operators in the 
effective Lagrangian obtained by turning on and off some coefficients (later indicated 
with M's) that regulate the weight of the operators. In fact, most of these coefficients 
turn out to be in principle free parameters (with the sole requirement for them to be 
smaller than the mass of the underlying theory) and this clearly enlarges the region of 
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the parameters space than can be spanned by the theory. Specific inflationary models 
would, on the other hand, put strong bounds on several of these coefficients. Besides 
this fact, the effective Lagrangian naturally captures extrinsic curvature-generated op- 
erators that are sometimes neglected in the literature. These terms should definitely be 
studied as, in fact, their contribution can be relevant [2^1125] and therefore they can also 
also significantly increase the dimensions of the parameters space of the theory. These 
are the specific operators this paper will mainly focus on. One more advantage that 
comes with employing the proposed setup is of calculational nature: in the so called 
decoupling regime (which implies to work in a specific energy range) the dynamics of 
the metric decouples from the one of the scalar that drives inflation thus rendering the 
Lagrangian itself and the higher order correlators much easier to handle and calculate. 
This mechanism is very reminiscent of what happens in standard quantum field theory 
and goes under the name of equivalence theorem. 

Interestingly, the effective approach proves itself useful already at quadratic order in the 
perturbations in that it automatically generates an equation of motion whose solution 
will encompass the classical wavefunction of many inflationary theories in the appropri- 
ate limits. We present here in detail such a solution. First, we give results at leading 
order in generalized slow-roll parameters for the power spectrum, its tilt n s and running. 
Secondly, we solve the equation of motion at next-to-leading order in slow-roll, then we 
specialize to the Ghost Inflation case and obtain tilt and running as well. 
The amplitude f^L has already been calculated in a number of papers [231 [26J [271 12H] 
within the effective field theory set up. Here in particular, we use the results of [21] 
to calculate the running of the non-Gaussianity (NG). From (almost) scale invariance 
arguments one expects the latter to be a weak effect but it is nevertheless quite impor- 
tant. In fact, one expects the scale invariance to be eventually broken and also wants 
to use the running [29]- [36] as a tool to distinguish between inflationary models which 
generate the same predictions, even at the level of the bispectrum shape-function. 
The paper is organized as follows: in Section 2 we solve at leading order the generalized 
equation of motion, we compare it with known solutions in the proper limits and com- 
ment on the wavefunction behaviour deep inside the horizon and in the crossing-horizon 
region; in Section 3 we give the general quadratic solution up to next-to-leading order in 
slow-roll, we then specialize to Ghost Inflation and calculate the tilt and the running of 
the power spectrum; in Section 4 we calculate the running of the f NL portion generated 
by an intriguing third-order interaction term introduced in [21], weighted by the M§ 
coefficient, and show that it can be significant; in the Conclusions we summarize our 
results and comment on further work; in the Appendix we give explicit expressions for 
several quantities whose form has been kept compact in the main text for the sake of 
simplicity. 
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2. Solution to the equation of motion for the scalar 

We here solve the equation of motion for the second-order effective Lagrangian derived 
in [23] (see also [2H [28] ) at leading order in slow-roll: 

C 2 = a\MlH{d^f + 2M 2 V + MfH^f + ^(d 2 n) 2 + ^(^) 2 ). (1) 
First though, a few comments on the above equation are in order. 

• The procedure according to which the Lagrangian was obtained is outlined in [23] 
and allows for a very general expression for inflation driven by a single scalar degree 
of freedom. In short, Eq. ([1]) is the most general second-order Lagrangian in unitary 
(co-moving) gauge provided the approximate symmetry of the underlying theory is 
such that only derivative terms of n appear in the action (see also [37] for an in- 
depth analysis of these issues). To first order, (which is all we need here) the scalar 
it is lineraly related to the dimensionless gauge invariant quantity £ by £ = —Hn. 

• In full generality the M coefficients above should be time dependent (we will deal 
with such a case in Section 3). However, if one is only interested in performing 
leading-order calculations, then, due to a generalized slow-roll approximation, one 
can safely consider them as constant. 

• Note that there is no trace of metric perturbations in Eq. (CQ). This is because we 
are working in the so called decoupling regime: for a sufficiently high energy range 
the dynamics of the scalar degree of freedom which drives inflation is decoupled 
from gravity. This is the so-called equivalence theorem at work. It suffices to say 
here that one can safely work with Eq. ([I]) assuming E > e 1 / 2 !!, where as usual 
e = H/H 2 , and E > M 2 /M Pl . 

• Note for example that for M 2 = = Mi,2,3 one re-obtains the usual quadratic 
Lagrangian with the speed of sound c 2 = 1 and the standard oscillatory solution to 
the equations of motion. On the other hand, switching on the M 2 term amounts 
to allowing c 2 s < 1 models, 1/c 2 = 1 - 2M 2 4 /(Mj v #), which have been proven to 
give larger non-Gaussianity [231 EH [391 SO, SI] • Working in the de Sitter limit and 
turning on M 2 ^ one rediscovers Ghost Inflation [12]. This very same procedure 
enables one to capture all the corresponding models at higher orders as well. The 
list of correspondences can be continued with K-inflation [131 H3] theories and 
others proving that the effective action approach naturally provides a more unifying 
perspective on inflationary models. 

• The action in Eq. ([TJ and, more importantly, its higher order counterparts, are 
generally written with large non-Gaussianities in mind. In the effective theory 
language, generalizing what has been done for DBI models, this corresponds 
to putting some constraints on the parameter space spanned by some time- 
independent (at leading order in slow-roll) coefficients in the action. A case in point 
is requiring a small speed of sound: this assumption, which can in general lead to 
to large NG, automatically translates into bounds on the values of the coefficients 
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driving quadratic operators in the theory Lagrangian. A speed of sound different 
than unity necessarily generates a different weight, for example in Fourier space, 
between time-like and space-like derivatives acting on the scalar. A straightforward 
generalization of this argument shows that, in writing each M, M coefficient driving 
a specific operator in the Lagrangian, there is a meaningful criterion for a selection 
as to which are the leading terms to write down. To illustrate how one goes about 
determining which contributions fit in the action, we now consider a simple example. 
Let's take terms that are multiplied by the Mf coefficient. In principle we should 
see in the above Lagrangian the contribution: 

-Ml(-H{d^f/a 2 + QHTx 2 ). (2) 

On the other hand, we will show below that, in comparing terms of the same 
perturbative order within the same M coefficient one just counts the number of 
space and time derivatives: the term with the highest number of space derivatives 
will be the dominating one. This is reminiscent of the fact that in DBI theories, 
when in the horizon crossing region, one can safely assume the following estimates 
to hold: 7r ~ Hit; Vtt ~ H/c s ir and so for c s C 1 space-like derivatives dominate. 
We will provide a generalization of this argument (see Eq. ([7]) below). Let us stress 
already at this stage though that, should one decide to include all these sub leading 
contributions in the action, the functional expression of the solution will not change, 
one merely redefines a couple of time-independent coefficients. This is due to the 
fact that the types of operator 7r 2 , (<9j7r) 2 , (<9 2 7r) 2 are already saturated at the level 
of Eq. (CD). 

We are now ready to tackle the equation of motion. After the usual change of variable, 
ir(k, t(r)) = a(r)u(k,r), the equation of motion can be written as: 
2 

u" - —u + a k 2 u + (3 k 4 r 2 u = 0, (3) 



r 2 

where «o, /?o are time independent (again, at leading order) dimensionless coefficients. 
This equation has been written in the context of tilted Ghost Inflation [26] and to our 
knowledge, it has not been solved analitically before Ref. [23], where the analytical 
solution has been briefly introduced and used for the computation of the three-point 
function. Here we discuss in much more details the properties of this solution. At this 
stage one can immediately recognize «o as the more common c 2 and Pq as the constant 
a 2 H 2 /M 2 first introduced in [12]. The complete expression for the coefficients is: 

= ~M 2 Pl H-MfH = (M 2 + M 2 )H 2 

° —MptH + 2M| ' 2(-M 2 l H + 2M$Y [) 

so that one reobtains the actual c 2 for Mi = 0. Note that one can simply look up the 
e.o.m. solution for DBI-like inflation if f3 = = Mi and Ghost Inflation in the de Sitter 
limit provided «o = 0. Let us pause here to comment on the possibility of a negative 
ao (see also [28]). Such a scenario would result in a region in the fc-space, whenever 

\a \k 2 > (3 k 4 r 2 - \, (5) 



for which the solution to the equation of motion will behave exponentially. Such a 
possibility raises a number of issues we address below. First of all, in order to keep 
control of the negative «o region of the parameters space of the theory in the ultraviolet, 
one requires that the (positive) f3ok 4 T 2 prevails over the «o contribution before k reaches 
the cutoff scale A. Considering that on the IR side, as we will show, the modes will 
eventually freeze outside the horizon, the case of a negative oiq should not in principle be 
disregarded. On the other hand, a lot of care should be exerted because an exponential 
phase of the modes for a sufficiently wide k region could generate values for higher order 
correlators that directly contradict available observational data. 

We could now proceed to solve the complete equation of motion. Equipped with just 
the equation fl3]), we can already make some educated guesses on the behaviour of the 
wavefunction. First of all, the typical oscillatory behaviour deep inside the horizon is 
to be expected in this more general well: both a^k 2 and (3ok 4 r 2 cause wave- 

like behaviour (see Fig 1 below) of the wavefunction, while the (— 2/r 2 ) contribution 
is negligible. This is important in that it tells us the main contribution to correlation 
functions will be coming, as usual, from the horizon- crossing region. Note here that, as 
far as (3o 7^ 0, the 'Ghost Inflation' term will eventually lead the oscillation if one goes 
deep enough inside the horizon. 

On the other hand, in the r — > limit, (—2/r 2 ) will be leading the dynamics and we 
expect to recover the usual, frozen modes. As is familiar from the DBI-like cases, it is 
convenient to introduce the notion of an effective horizon, placing it where the oscillatory 
behaviour stops being dominant. In formulas: 

a k 2 + (3 k 4 r 2 = A => n = 2 (6) 



T 



k\ oiq + V"o + 8$ 



For (3q = 0, ao ~ 1 one recovers k 2 r 2 ~ 1 at the horizon. 

At this stage we can perform a consistency check and show how one can generalize the 
argument, initially borrowed from DBI-like inflationary models, that in comparing terms 
at the same order in perturbations and with the same overall number of derivatives, 
the ones with the most space derivatives are dominating in the c s C 1 limit. The 
generalization of this argument consists in restricting the parameters space to the 
ao < 1 and /?o <^ 1 region. Consider Eq. ((3]) in Fourier space; in full generality 
one expects V7r ~ kn and n ~ Hit so what needs to be done is relate k with H at the 
horizon. Using equation ([6]) and r ~ —l/(aH) one obtains 

k = V2H 



«o + v«o + 8 A) 

Since the main contributions to correlators comes from the horizon- crossing region, 
this shows that, for (a 0) Po) < 1 we can still identify leading terms in the Lagrangian 
according to the standard procedure. 



6 



2.1. Wavefunction 



Let us verify all this quantitatively. The solution to Eq. (j3]), being of second order, will 
come with two ^-dependent integration constants. We have determined their values by 
requiring to re-obtain the known DBI and Ghost solutions in the corresponding limits. 
The general wavefunction reads: 



u k (r) 



2V^ 



1 

4 



IQ>0 



Ci(fc) 



+ 



2V4, 



c 



1 , ia 3 f— 2 2 

4 + 4^'-2'- Z V^ T 



C 2 (fc) 



(8) 



Where Q stands for the confluent hypergeometric function and C is the generalized 
Laguerre polynomial. We verified that, properly adjusting the integration constants 
according to 



Ci(fc) 



-\ -3/4 

) r 


5 iao 




L 4 



fc-3/2 



M|,e H + 2M$ 2V4 r 



3/2 



4 iao+V^o 



C 2 (fc) = 0, 



(9) 



one obtains, in the appropriate limits [45] , the wavefunctions of standard inflation and 
Ghost Inflation [42]. We can now write our solution: 



7Tfc(T) 



4 v po 



VM2eF 2 + 2M2 4 y2 7o 3/4 r(| + 



3/4^^5 j qq_ 

4 Aao-Uyjfc' 



(10) 



where 7 = ao + a/A) and is the Euler gamma function. 

We see in particular that for ao = 0, Eq. ( fTUj) immediately reduces analitically to the 
Ghost Inflation wavefunction 7r fc (r) = (H (-r) 3/2 / v / 2M|) ^/f Ul^^^fW^r 2 ) with 
being the Hankel function of the first kind. On the other hand, one can easily see nu- 
merically that the DBI solution is recovered in the /? — ?• limit. 

To give some intuition on the behaviour of the general, interpolating wavefunction, we 
plot it in several (ao, 0o) configurations. For overall consistency in the comparisons, in 
all the following pictures we have chosen points in the (ao, /3o)-pl ane so that the horizon 
crossing always lies at the same point, numerically r* = — a/2, and we have plotted the 
wavefunction from well inside the horizon (r = —10 r*) up to r = 0. 
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Figure 1. On the left, the general wavefunction in the DBI-likc, /3q —> 0, limit; on the 
right the DBI-like solution itself. We plot the real part and find perfect agreement, 
same holds for the imaginary part. To produce the plot the parameters have been set 
to: ao = 0.1, Pq — > 0, k = 1, H = 1 and the Planck mass-dependent normalization has 
been neglected. The corresponding plot can be omitted for the Ghost limit since in 
that case we recover the Ghost solution analytically. 




Figure 2. On the left, the Ghost Inflation wavefunction (ao = 0). On the right the 
general interpolating solution calculated for ao = 1/2 , (3q = 1/4. 



From these plots we learn several things. First, as argued before, there is a common 
oscillatory behaviour once inside the horizon. The frequency of these oscillations is 
more pronounced for the Ghost solution when deeper inside the horizon. In the general 
solution the frequency varies according to the "relative weight" of the Ghost component, 
Po, and the DBI-like one, ao- 



2.2. Power spectrum 

We now turn to the expression of the power spectrum 



P, = ^\n(k,r^0)\ 2 - 



2tt 2 1 v ' n ' 16ir(Mj,eH 2 + 2M*)(a + v %) 3/2 \T(i + 7-^ 



(11) 

Clearly there is no time dependence in the above result, the modes freeze outside the 
horizon. Note here that, to reproduce standard results, we should re-introduce the speed 
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of sound c s which, in the language we are using, is related to M 2 via: 
1 2M 4 



^ = 1 2 



c 2 s M 2 P H 



Upon switching to the gauge invariant quantity related to 7r by £ = — Hn [23] (see also 
the appendix of [27] for the relation at second order), and reintroducing the proper units 
with Planck mass we get: 

p = {a ° + ^)~ 3/2g4 no) 

C 1 MM^ + 2M|)|r(| + i -^= ) p- [ } 

Again, one could easily check that the above result analytically covers the power 
spectrum of DBI-like and Ghost Inflation. 

As we mentioned in the last comments to equation (pQ), even when including sub leading 
terms in the quadratic action, the functional dependence of our wavefunction does not 
change, only the definition of a , (3 does. Since we are now going to set bounds on 
operators coefficients , we want to be as precise as possible and will therefore extend the 
definition of the two parameters to cover the subleading terms as well. We now have: 

a = -MhH - M?H/2 = M 2 H 2 /2 

-M 2 Pl H + 2M 2 4 - 3M?H ' -Mlfi + 2M 2 4 - 3M?H ' 

In obtaining Eq. ( f!3|) . we took into account the fact that the Mi, M 2 , M3- driven terms 
multiply operators of the type tt 2 , (diii) 2 as well. We also choose to replace M 2 and M 3 
with a linear combination of the two masses: we set Mf = -3M| and M 2 = Mf + M|, 
see also [23] ■ This procedure allows one to put to zero all the subleading operators 
tuned by M 2 ,M3 and makes the correspondence between inflationary models and the 
switching of the M, M parameters absolutely sharp f. 

An immediate simplification is that now, upon requiring M = -v^> f3 = 0, one goes into 
DBI inflation, exactly. Similarly, now de-Sitter limit and Mi = give Ghost Inflation 
with a = 0. The power spectrum looks very similar to the one in Eq. (fT2f) 

P = (<*«> + v^r 372 ^ 4 (u) 

C ySixiMp^H 2 + 2M 2 — 3MiH) + ^ ^ l2 ' 



4a -4iWI3o 



only the definitions of the «o, Po parameters and the normalization constant have 
slightly changed. Expressing all the parameters in the spectrum in terms of the M, M 
coefficients, we count the degrees of freedom as being five, associated to H, e, M 2 , M , M\. 
The first three are the same that appear also in the standard case as H, e, c s ; using c s 
instead of M 2 is just a matter of dictionary. In the Ghost Inflation case, the quantity 
Mo replaces the speed of sound and we are again back to three parameters. In the 
most general case one has to keep both M and Mi as well. Bounds can be put on the 
values of these five parameters by employing the expected value for the power spectrum, 
Pfi ~ 10 -10 and its tilt. Let us also mention here that further mild inequalities must 

I The reader might worry that one degree of freedom is lost. However, the two coefficients multiply 
basically the same interaction terms in the action up to fourth order. 
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be satisfied by the a , (3 parameters in order to keep the generalized speed of sound 
small, see [23]. In Sec. 2.3 below we present a calculation for the tilt and running 
of the power spectrum. These quantities are essentially obtained from considering the 
time dependence of the M coefficients in Eq. flH}. On the other hand, as has been 
specified above, part of the procedure that led to Eq. (fill) has been to disregard the 
time dependence of said coefficients §. Restoring it at a later step, as we do below, 
is standard accepted procedure because generally only at this stage the effect of time- 
dependence becomes important. A calculation that does without this assumption is 
presented in Sec. 3. 



2.3. Tilt and running 

Below we employ some simplification in order to present our result for the tilt of the 
power spectrum in a way that resembles as closely as possible the typical expression for 
n s — 1. Indeed, the spectrum dependence on the Euler Y function in Eq. (fl~4")) is not to 
be found in e.g. DBI, Ghost Inflation etc. For simplicity, we choose not to write here 
the explicit dependence of the Euler function on the M coefficients and leave it implicit; 
we report the full explicit dependence in the Appendix. We apply the following formula, 



d /dink 
n. — 1 = — — - In Pi 



-i 



d In k \ dt 



1 dP h 



1 dPv 



*=*• HPt dt 



*=*7 Pk dt 

on the power spectrum, where t* is the time at horizon crossing. The time dependence 
of the M, M coefficients is taken into account and the time dependence of the Euler 
function is dealt with as one would do with a generic function T(t). One obtains: 

f 7H 2 M 2 e + 8 (2M 2 4 + 3M! 4 ) 

Us ~ 1 ~ ~HT ~ E X 2 {H 2 M 2 P e + 2M 2 4 + ZM^) + 
8H 4 M P e 2 - 16H 2 M 2 e (2M 2 4 + 3M^) + H 2 M 2 e (-UH 2 M 2 e + 8M 2 4 + 15M! 4 ) 

6 X 



2 (H 2 M 2 P e + 2M 2 4 + 3Mx 4 ) \ 2H 2 M 2 P e + M x 4 + M 2 ] /2H 2 Mpe + 4M 2 4 + 6M^ 

(QH 2 Mpe + 24M 2 4 + 33Mi 4 ) H 2 Mpe 



4 (H 2 Mpe + 2M 2 4 + 3Mi 4 ) [2H 2 Mpe + M x 4 + M 2 J2H 2 Mpe + 4M 2 4 + 6M X 4 



— 7] X 



H 2 M 2 e M 2 2M 2 4 

~ V X 4 (H 2 Mpe + 2M 2 4 + 3Mx 4 ) ~ HM 2 X (H 2 Mpe + 2M 2 4 + 3M X 4 ) + 

M 2 (QH 2 Mpe + 3Mi 4 ) 2M 2 4 

H x 

HM. 



2 (H 2 Mpe + 2M 2 4 + 3Mx 4 ) \ 2H 2 M 2 P e + M x 4 + M 2 ] /2H 2 Mpe + 4M 2 4 + 6M X 4 
Mi (+3M X 4 - 4H 2 M 2 e + 4M 2 4 ) 3M X 4 

X 



HMl (H 2 Mpe + 2M 2 4 + 3Mi 4 ) \ 2H 2 M 2 P e + Mi 4 + M 2 J2H 2 Mpe + 4M 2 4 + 6M : 4 



§ This is because, at the level of the action, considering the time dependence of these coefficients would 
automatically translate into going at next-to-leading order in slow-roll; for a leading-order calculation 
it is therefore sufficient to consider their values at the horizon. 
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M 6M 2 



T x 7 — 

#M / 2Mo 2 + ^2 (2H 2 M 2 e + M x 4 ) 



_ff 2 Af£ e+2M 2 4 +3Mi 4 



M x 3M1 4 



FMi (H 2 M 2 P e + 2M 2 4 + 3Mi 4 , 
where all the quantities are to be intended as calculated at horizon crossing. 
For convenience, we have factored out the usual parameters: e, 77, s (the latter is written 
in M 2 language; we give the dictionary in Eq. (ITS]) below) and their generalization: 

M Mi f 

e °-HW ] 61 ~ ~HM\ ' er "OT- (17) 

£2 ^^ = ^-e+ « (18) 
#Af 2 2 2( C 2 - 1) v ; 

The variable c s is the one defined in Sec. 2.1 and 77 = e/(He). For simplicity, we have 
defined the variables M 4 = M$H 2 and M x 4 = Mfif . 

Let us briefly comment on the above results. First note that, as expected, the factor that 
each of these generalized slow- roll parameter multiplies is of order unity or smaller. A 
quick consistency check consists in specializing the formula above to known inflationary 
models, for example requiring Mi = = Mo gives back the usual result [27] for DBI-like 
models: 

n s - 1 = -2e - 77 - s (19) 

Consider now the more general case with (Mi = , H ^ , Mo 7^ ), which comprises 
DBI-like theories and Ghost Inflation models as limiting cases. In such a scenario there 
are some mild bounds to be required on Mq , MpeH 2 . First of all, since we are interested 
in the «o 1 , (3q 1 region of the parameters space this requires: 

M 2 P eH 2 < M 2 4 ; M 4 < M 2 4 . (20) 

On the other hand, the M -driven slow-roll parameters is eo = M /(HM ) and therefore 
the above inequalities do not put upper bounds on this slow-roll parameter. Much 
like we will see in the next section for the running of the bispectrum amplitude, one 
must instead be careful to account for the fact that too large a value for Mo could 
give a contribution to the bispectrum amplitude that must be excluded. Indeed, once 
expanded, the Mo-proportional contribution to the quadratic action for the scalar reads: 



[ d 4 x^[. + M 4 (t + tt)^-} ~ / d A x^[. + M 4 (t)^^ + 4M3(t)M 7r 

(21) 

with the second term on the RHS action clearly contributing to the three-point function 
of the scalar. This M -generated contribution must be weighted against the third order 
interactions generated by M itself, but also by new M coefficients that first appear at 
third order in the action (see [21] and the analysis in Sec. ^). A similar analysis applies 
for e 1 if we let M 1 ^0. 



a 4 
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We can conclude that the mild bounds on these new, generalized slow-roll parameters 
come from the obvious fact that we are doing a slow-roll expansion, from the value of 
the power spectrum itself and from the requirement that they are not so large as to 
produce too large a value for f NL . 

In order to obtain the running of the power spectrum we proceed as below: 
dn s 1 dn s 

as = dh^k ~ H~df- ( } 

We give below a compact results: 

r 2 e ® ^©2 £202 V® v yQ v eiQi eiQi epQo e 9 

a s - H2T e r H H H H H H H H H H ' 

(23) 

and point the reader to the Appendix for a more explicit expression of the coefficients 
functions 0, 2 , 0,,, ©i, 0o- 

3. Next-to-leading order 

The discussion presented so far is based on a generalized slow-roll approximation at 
leading order. In particular, the M, M coefficients driving the various operators in the 
Lagrangian are assumed to be time independent. This assumption propagates into the 
equation of motion for the scalar ir, the classical solution itself and the power spectrum. 
To make up for this approximation when calculating the the tilt of the spectrum, one 
restores the time dependence of the coefficients at the level of the power spectrum. 
A more systematic approach consists in accounting for the time dependence of the 
coefficients already at the Lagrangian level, this is done by taking the generalized slow- 
roll approximation to next order. Schematically one has: 

oc J d 4 xy/=g[ - Mp(3H 2 (t + tt) + H(t + tt)) + M(t + tt) x (quadratic)] . (24) 

The first term does not appear in the Lagrangian in Eq. ([1]) because, at leading order in 
slow roll, it is not quadratic in fluctuations, it does in fact, contribute to the background. 
We see that accounting for the 'V in the time dependence of M in the last term in 
Eq. (I24p would result in a third order operator. This must be considered when studying 
interactions but it is not what we want to analyze here, the wavefunction comes from 
the quadratic Lagrangian. On the other hand, the 'V in the first term of the action has 
to be accounted for; doing so results in just one additional contribution to the action 
and it turns out to be proportional to e 2 . This is all consistent with the fact that, at 
leading order, the action is instead proportional to e. We now have: 



Co = a 



M 2 P H(t)(d^y + 2M 2 4 (t)7r 2 - Mf(t)H(t)(3r 



2a 2 ' 

t) Vl . (25) 
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We proceed to write down the equation of motion as in the leading order case, obtaining: 



o k + a (l + s a )k 2 a k + (3 (1 + sp) 



a 2 H 2 



<7k 




M 2 H(t) 



2M 2 4 - M 2 P H - MfH , 



(26) 



where the following definitions have been employed: 



H(Xq ^ 



sp 



MpH + 2M 2 4 - 3M?H) 



7T 



a 
1' 



(27) 

In order to solve the equation of motion one needs to calculate /"// explicitly; the result 
is given in the Appendix. A compact expression for the equation of motion is given by: 



a k + a k a k + (3 k r a k = (1 + x )a k 



{2i 



where Xq is a linear combination of slow roll parameters and &o, (3 represent a slight 
redefinition of the initial parameters. For explicit expression we refer once again the 
reader to the Appendix. Equipped with Eq. ( 1281) . one uses the Bunch-Davies vacuum 
condition to write down the solution: 



a{x,k) = Ci(fc)- 



-ia 



/3o+2/3 +/3oV 9+8a;o i 
' 2 



4/?o 



2 + V9 + 8x ) , -i\/(3 k 2 x 2 



2 -i(2+ v / 9+8^) e -I 



,2^(2+v / 9+8^) 



(29) 



(X*)' 1 ^ ' \JX 

where Q is the usual hypergeometric function, x = — r, and the wavefunction must be 
expanded to first order in xq = 0. In the first section, we gave the exact expression 
for the solution above at leading order and went on to calculate the resulting power 
spectrum, its tilt and running. The calculation at next-to-leading order has already 
been performed for DBI-like theories of inflation, using the same formalism employed 
here, in [27]. Here instead, we choose to calculate the next-to- leading order Ghost 
Inflation solution obtaining also the tilt of the power spectrum and the running. The 
procedure is a standard one, so we briefly sketch it. The two ^-dependent constant of 
the Ghost equation of motion 

a'i + j3 k\ 2 a k = — (1 + x%)a k , (30) 

are reduced to one by imposing the correct leading order limit on the wavefunction. 
The remaining constant is obtained by requiring the proper normalization, that is by 
imposing the following commutation relations to hold: 



7T (x 



P 



i5\x-y); [a % ,al\ = (2n) 3 5 3 (k - p) , (31) 

where P is the momentum conjugate of the scalar 7r and the creation and annihilation 
operators a, a' are the usual operators in the free field expansion for the quantized field 
7r. Proceeding as prescribed above, one obtains: 



71 



(k, 



H 



2tt(1 + spfl\-rfl 2 H<V [|V9T8^, \k 2 ^ (1 + s P Y 



3 1 + ? 



X { 



7/3-f + log2) r l + I v /9 + 8 



ix 



(32) 
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The wavefunction above is the correct Ghost Inflation wavefunction up to next-to- 
leading order in generalized slow-roll parameters; it gives back the leading order solution 
and respects the proper Bunch-Davies vacuum requirement. From Eq. ( 1321) one can 
readily calculate the tilt of the spectrum: one simply considers the leading behaviour of 
the wavefunction as k goes to zero. In fact, upon expanding for small k one finds that 
our solution goes like: 

"3 x 1 /— ( sf 



A ' 3 ' 

from which we obtain that 



k 



2_2 



3 2a:n 

k-i-^ 1 , (33) 



n s - 1 = -4/3x - (34) 

When specialized to Ghost Inflation, the value of the xq parameter (which always 
constists of a linear combination of the generalized slow roll parameters) is given by: 

G n M 2 

X ° = (35) 

and the running amounts to simply 



din A; 

4. Running of /jvl 



Ax /3H. (36) 



In the quest for properties that help in removing degeneracies among the many 
inflationary models one generally considers another observable beyond the power 
spectrum and its running, i.e. the analysis of non-Gaussianities. Starting with the 
bispectrum, one can study its amplitude, shape and running. In the same spirit of the 
analysis we performed for the power spectrum, we now want to estimate the value for 
the running of the bispectrum amplitude, /jvl- In performing the calculation for P^, 
n s — 1 and then a s , we used the fact that the coefficients driving quadratic operators in 
the Lagrangian are nearly constant, up to slow-roll corrections. At first approximation, 
one writes down the power spectrum as a function of these parameters calculated at the 
horizon. Only when calculating the tilt of the spectrum and its running one does consider 
the time dependence on the M's, thus obtaining Eq. (116j) . (123j) . Similarly here, we will 
employ the results on the bispectrum amplitude contributions generated by independent 
interaction terms as given in [21] (see [17]- [50] for the in-in formalism). In particular 
we are going to focus on the running of /jvl generated by a third-order interaction 
term whose bispectrum shape-function peaks is an uncommon flat configuration (see 
also [25]). In fact, in [21] there are several independent terms that generate at least two 
qualitatively different flat shape-functions; we choose here to concentrate on the analysis 
of just one term as the same considerations and conclusions can be straighforwardly 
adapted to all of them. 

The analysis for the running of /jvl has been done for several inflationary mechanisms 
such as DBI inflation, and others. We are going to extend this type of study to a 
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third order interaction term driven by a nearly constant coefficient, M 6 (t), and which 
is generated by an extrinsic curvature contribution: | 

J d 3 xdty/^[... + 6K } - E(^'^) 2 + •••] • (37) 

hi 

In [2lj we calculated the corresponding contribution to f^L in six different points of the 
(a , (3 ) plane, we report them in Table 1. 



Benchmarks 


1 


2 


3 


4 


5 


6 




1(T 2 





0.5 • 10" 2 


2 • 10^ 7 


10~ 4 


10~ 6 


00 





0.5- 10" 4 


0.25 • 10^ 4 


5- 10" 5 









Table 1 : configuration 1 describes pure DBI-likc theories, pure ghost corresponds to configuration 2. In 3,4 a more 
general model is considered while in the last two configurations one aims at considering the cases characterized by very 
small generalized speed of sound, ^/oq. 



Corresponding to each one of the configurations of Table 1, we report in Table 2 below 
the value of /nl- 



benchmarks 


1 


2 


3 


4 


5 


6 


fM 6 
J NL 


10 4 76 


4-10 3 76 


5-10 3 7 6 


4-10 3 76 


10 8 7e 


10 12 7e 



Table 2: the value of the bispectrum amplitude f^fl m the six different configurations. 

In Table 2 we introduced the dimensionless quantity 76, which is given by 

7 6 = (M 2 H 2 )/(M Pl eH 2 + 2M 4 - 3M 3 if). (3 

Computing the running of f^L amounts then to calculating 

d\n\f NL (k)\ 1 df NL 



NG dink HfNL dt 

For the specific case at hand, we have: 



(39) 



1 d 



1 d r 



fM 6 _ 

HfS[dt lNL HfSldt 



4M 2 4 + 2Mi 



t=t* 



(H 2 M 2 P e + 2M 2 4 + 3M X 4 ; 



H 2 M 2 e 



8M0 



'' H 2 M 2 P e + 2M 2 4 + 3Mi 4 62 (H 2 Mpe + 2M 2 4 + 3M X 4 ) ^ (H 2 Mpe + 2M 2 4 + 3Mi 4 ) ' 

(40) 

In the above expression N stands for a numerical factor which is dependent on ao and 
f3 , specifically it goes like (a + A))~ n with n small positive integer. The generalized 
slow-roll parameters are defined as usual, with the only new ones being 66, e^: 

M 6 



+ 2e 6 + e^v + 
4Mi 4 



e 6 



e N = 0(1) x (e, 77,62,61,60) 



(41) 



The contribution of Me to the action at third order is proportional to — 7\/|(g 00 + 1)8K"8K£ 
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where K^ v is the extrinsic curvature tensor 



The last expression in Eq. (141 j) means that e^v can be expressed as a linear combination 
of the generalized slow roll parameters introduced before multiplied at most by an order 
unity constant. Overall, we can then conclude that the running of f N [ is of the order of 
the generalized slow roll parameters (or a linear combination thereof). Now, all these 
parameters, except for e^, can be expressed as functions of M, M coefficients that first 
appear in the quadratic Lagrangian of the theory. As such, four of these parameters 
could in principle be expressed as a function of the observables one usually uses, namely 
the scalar and tensor spectral indices and their running. Not so for 66, on which, in 
principle, we enjoy more freedom as it drives terms in the action that are at least cubic. 
It could indeed be that the running of f^f[ is dominated by this contribution and could 
therefore be larger than one finds in some general single-field slow-roll models where 
the extrinsic curvature-generated interaction terms are not accounted for. The same 
considerations apply to other interaction terms such as the Mg-driven one in [21]. On 
the other hand, special care must be exerted so as to make sure that requiring e 6 to 
be the leading generalized slow-roll parameter in Eq. f PTOl) does not spoil the possibility 
to have the contribution in Eq. fl37j) dominate the overall bispectrum amplitude, which 
is what made this type of contribution interesting in the first place. Indeed, there are 
two ways of making e 6 the leading parameter, a large M e and a small M 6 . Pushing the 
latter option too far the interesting and possibly leading bispectrum and trispectrum 
amplitudes [21] would become subdominant and this would make the corresponding 
flat shape-functions a mere curiosity. But also the former option has to be discussed 
and this is clear from the expansion of Eq. ( 1371) resulting from considering the Mq time 
dependence: 



/ 



, . r M 6 (t) 2 7T x9 2M 6 (t)M 6 (t) n v^.„ ,o i 

O CL ■ ■ O CL ■ ■ 



Indeed the quartic interaction term that appears is proportional to M s and if the latter 
is too big it could give rise to too large a contribution to the power spectrum at one 
loop and would have to be ruled out. Below we give a number of inequalities that 
the quantities M 6 , M 6 need to satisfy in order not to spoil the appealing bispectrum 
and trispectrum features outlined above. It turns out they are not too restrictive and 
that the running of fjtfl can be safely ruled by €q. We first write down the inequalities 
that stem from requiring that the Mg-proportional quartic interaction term is not the 
leading interaction in the fourth order Lagrangian (this would have consequences on the 
loop corrections to the power spectrum as well) as compared to the usual single-field 
interactions (for the complete action at fourth-order see Ref. [25] ) 

K >e«»e; {a °:f?^>e 6 »e; ^^>e 6 »e. (43) 



M|# 2 - ^ " ' M2 H 2 M 2 H 2 

Notice that only one of these inequalities need be satisfied. The e$ ^> e part of 
the inequalities above ensures that indeed e$ is the dominating generalized slow-roll 
parameter. One has to keep in mind here that large non-Gaussinities are generated by 
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requiring the generalized speed of sound, ~ « + a/^, to be much smaller than unity and 
so the first of these inequalities seems somewhat less stringent than the others though 
there is no requirement on the M's to be all of the same order. 

As anticipated a small M 6 can in principle lead to a subleading contribution to the 
bispectrum signal thus rendering the corresponding flat shape-function less interesting. 
Borrowing the third order action of [241 |25] and employing the estimates on the 
wavefunction we showed to hold at horizon crossing, one is able to derive the inequalities 
below: 

> 1 ; ,,4, 6 , 7^ > 1 • (44) 



Mi(a + V£o) ' M 3 > + vOT 
Again, it is sufficient that only one inequality in Eq. ( 143]) and the corresponding one 
in Eq. fj44|) are satisfied. The coefficient M 4 is not found in the equation above as it 
first appears in theory Lagrangian at fourth order, so it is not involved in the tree level 
bispectrum calculations. The further freedom on M 4 (and other coefficients) that results 
from this simple fact can be used to study models of inflation which present a relatively 
small bispectrum together with a larger trispectrum signal [251 US]- The inequalities 
given in Eq. ( 1431) and ( )44|) above are indeed compatible for values of («o + a/A)) smaller 
than unity, an assumption which is generally made when looking for models that can 
produce large NG, as shown in [241 [25] . We see then that there is a whole, large window 
of values for the coefficient M e that would allow for a running of f^L dominated by 
€q. By looking at Eq. ( 140]) one realizes the running can indeed be large since it was 
shown in [24] that itself can give the leading contribution to the total bispectrum 
amplitude. This effect is generated purely by third order terms and specifically by 
extrinsic curvature-generated interaction terms driven by coefficients which first appear 
in the Lagrangian at third order. Adding the results of this section to the analysis 
of [241 125] one can safely say that, with respect to all observables one is ultimately 
interested in, the analysis of these curvature terms has shown they can have leading 
effects on all quantities and must therefore always be included in a thorough analysis of 
non-Gaussianities. 



5. Conclusions 



In this work we employed the tools of effective field theory to analyze a very general 
theory of inflation driven by a single scalar degree of freedom. We have done so at the 
level of the quadratic and higher-order Lagrangian. The solution to the equation of 
motion for the scalar has been analyzed in detail. The freezing of the /c-modes outside 
the effective horizon is shown to be preserved in this setup as well. In Section 2 we gave 
the resulting power spectrum at leading order as a function of all possible parameters. 
We learned that, in order to consider a fully general theory, one has to raise the number 
of generalized slow-roll parameters to five, as opposed to the usual three (H, c s , e). 
Naturally, this fact allows for more freedom and we stress here that the two additional 
degrees of freedom (eo and e\, Eq. (1171) ) are generated by extrinsic curvature terms in 
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the quadratic Lagrangian that can be significant in specific cases such as near de Sitter 
limit, or for small values of the generalized speed of sound or a combination thereof. 
From the spectrum, we proceeded to obtained its tilt and running. In Sec. 2.3 we 
also discussed the (mild) bounds on the generalized slow roll parameters. In Sec. 3 
we described a next-to-leading order calculation for the tilt and running of the power 
spectrum, specialized to the Ghost Inflation case. At third order, we studied the running 
of the bispectrum amplitude f^L contributed by a curvature-generated interaction term 
driven by Mq. From previous investigations [24"l [25] . it was known that this term is 
able to generate important or even leading non-Gaussianities and presents distinctive 
features in the form of the bispectrum and trispectrum shape-function. In the last 
section we showed that, upon imposing suitable and mild bounds on the coefficients 
driving the interactions, the running of Jnl is dominated by this term without spoiling 
any of the aforementioned interesting features on the three and four-point function and 
without affecting the leading value of the power spectrum. In particular, we showed 
that the running ung can be such that ting 3> 0(e, rj). 

All the above results clearly point to the importance of considering extrinsic curvature 
generated terms in the quest for predictions on important observable quantities. 
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Appendix 

Explicit expression for the slow-roll parameter er 

We give here an explicit expression for the time dependence of the slow-roll parameters 
which we called er- This quantity if first written in terms of the M, M parameters: 



T[t] = 7 



H 2 Mp€ + iM x 4 



(H 2 M 2 P e + 2M 2 4 + 3Mi 4 ) 



4(_ff 2 Mj,e+iMi 4 ) 



Mi) 4 



X 



7 



+ 



H 2 M 2 P e 



(H 2 M 2 



9M4 + ( < H2M ^ +l 2^ 4 ) 

ZM 2 +6Mi ) I #2^+2^4+3^4 



(45) 



Mq 4 



2<H 2 Mle+2M 2 4 +3M 1 4 ) 
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where 7 ft] is the Euler function. From here one calculates the quantity er = — j^p , 
obtaining: 



HT 



\ 



PolyT 



\ 



' 4 



/ 8* (H 2 M 2 P e + 2M 2 4 + 3Mi 4 ) 



Mo 4 



2H 2 M 2 P e + Mi 4 



% [2H 2 M 2 P e + Mi 4 ) - 2 [H' z M'pe + 2M 2 4 + 3M : 



r2 it /r 2 . 



Mn 4 



( 8t (H 2 Mpe + 2M 2 4 + 3M X 4 ) ./^ 



2# 2 M|e + 4M 2 4 + 6Mi 4 i 
. -r 



Mq 4 



i£+4M 2 4 +6M 1 4 



2H 2 M 2 P e + Mi 4 



x 



i ( 2H 2 M 2 P e + MA + 2 (H 2 M 2 P e + 2M 2 4 + 3M X 4 



Mn 4 



2H 2 M 2 P e + 4M 2 4 + 6M1 4 



2H i M%e 2 



M n 4 



pC dt I 2# 2 M 2 e + 4M 2 4 + 6M! 4 ) + 4jff M P eM 2 ^ ( 2 #2 M 2 e + 4M 2 4 + 6M! 4 



r/ 



M 4 



+7ff 2 M|,eM 1 4 — 



+3M! 8 4- 



dt \2H 2 M 2 P e + 4M 2 4 + 6M! 4 
M 4 N 



2M 2 4 Mi 4 — 



d* \2H 2 M 2 P e + 4M 2 4 + 6M X 4 



dt V 2H 2 M 2 e + 4M 2 4 + 6M! 4 



lQHM 2 eM 2 4 H 



Mn 4 



^OHMpeM^H 



M 4 



m 2 M 2 P M 2 ^e 



2H 2 M 2 e + 4M 2 4 + 6M X 4 



2H 2 M 2 P e + 4M 2 4 + 6M X 4 



M n 4 



2H 2 Mle + 4M 2 4 + 6M1 4 



-10H 2 M 2 



M 



pMi 4 e ^ — 7 + 32# 2 M 2 eM 2 3 M 2 — = — 

p 1 2H 2 M 2 P e + 4M 2 4 + 6Mx 4 p 2 2 2H 2 M 2 P e + 4M 2 4 + 6M X 4 



M n 4 



+ 16M 2 3M 1 4 M 2 2g2M , e + 4M24 + 6Mi4 + 40H"^fm e+ "^ ^ 



Mn 4 



2 a #2 



M n 4 



Mn 4 



M 3 (4 (H 2 M 2 P e + 2M 2 4 + 3Mi 4 ) M - M (#M| foeH + Hi) + 8M 2 3 M 2 + 12Mi 3 M 



2 (H 2 Mpe + 2M 2 4 + 3M X 4 ) 2 



2M n 4 



2H 2 M 2 P e + 4M 2 4 + 6M1 4 
H 2 Mpe + 2M 2 4 + 3Mi 4 ^ 



(1 4 ) 2 x 



({2H 2 M 2 P e + Mi 4 ) 4 + 8 (2H 2 M 2 P e + Mi 

+ 16 (H 2 M 2 P e + 2M 2 4 + 3Mi 4 ] 



Mn 4 



2H 2 Mle + 4M 2 4 + 6M X 4 



M n 4 



2H 2 M 2 P e + 4M 2 4 + 6M X 4 



Explicit expression for the G functions of Eq. 



e 



7H 2 M 2 e + 8 (2M 2 4 + 3M X 4 ) 
2 (H 2 Mpe + 2M 2 4 + 3Mx 4 ) 



+ 
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&H A M A P e 2 - 16H 2 Mpe (2M 2 4 + 3Mi 4 ) + H 2 M 2 P e (-UH 2 Mpe + 8M 2 4 + 15Mr 



2 (H 2 Mpe + 2M 2 4 + 3M! 4 ) \2H 2 M 2 e + M x 4 + M 2 ^j2H 2 M 2 e + 4M 2 4 + 6M X 

" A (# 2 Mj,e + 2M 2 4 + 3M X 4 ) 

(QH 2 Mpe + 24M 2 4 + 33M! 4 ) # 2 M£e 



4 (H 2 Mpe + 2M 2 4 + 3Mi 4 ) [2H 2 M 2 P e + Mi 4 + M 2 ^/2H 2 Mpe + 4M 2 4 + 6M X 4 

e 2M/ | 

2 (# 2 M 2 e + 2M 2 4 + 3Mi 4 ) 

(<oH 2 M 2 e + 3Mx 4 ) 2M 2 4 



(H 2 M 2 p e + 2M 2 4 + 3M! 4 ) 2H 2 M 2 P e + Mi 4 + M 2 J2H 2 Mpe + 4M 2 4 + 6Mx 



i 



3M i 4 | 

1 (H 2 Mpe + 2M 2 4 + 3Mi 4 ) 



;+3Mi 4 - AH 2 Mpe + 4M 2 4 ) 3Mi 4 



(H 2 M 2 e + 2M 2 4 + 3Mx 4 ) 2H 2 M 2 e + M x 4 + M 2 J2H 2 M 2 e + 4M 2 4 + 6M! 4 



eo = r (47) 

(2AV + V2(2HV p e + AV) 7=^=^=^) 

Explicit expression for variables in Eq. (2b]). (2^) . 

f_ = 2a 2 H 2 L_ i , 6M 2 3 M 2 , 9Mi 3 Mi 



/ V 2 H ( 2M 2 4 + 3M i 4 - MpH) H (-2Af 2 4 - 3Mi 4 + MpH 

— mi. ^ , . (4f 

AH ( -2M 2 4 - 3M X 4 + MpH 



a = a (l + ^-) ; A, = A>(1 + ^) ■ (49) 
iiao Hpo 

In going from Eq. (1261) to Eq. (1281) it was safely assumed that a(r) ~ — Hr ^_ e y which 
in turn means that the cmplete and most general expression for the parameter x is, to 
first order in generalized slow-roll parameters, given by: 

e 6M 2 3 M 2 9Mi 3 Mi 
= I — — H 7 ; ; t^t + 



2 H (2M 2 4 - 3Mi 4 - M 2 p H) H (-2M 2 4 + 3Mx 4 + MpH 



AH (-2M 2 4 + 3M X 4 + M 2 pH) 2 2M 2 — MpH + 3M 4 
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